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Abstract. We study the intersection points of a fixed planar curve F with the nodal 
set of a translationally invariant and isotropic Gaussian random field ^(r) and the 
zeros of its normal derivative across the curve. The intersection points form a discrete 
random process which is the object of this study. The field probability distribution 
function is completely specified by the correlation G{\r — r'\) = (^(r)^(r')). Given 
an arbitrary G{\r — r'|), we compute the two point correlation function of the point 
process on the line, and derive other statistical measures (repulsion, rigidity) which 
characterize the short and long range correlations of the intersection points. We use 
these statistical measures to quantitatively characterize the complex patterns displayed 
by various kinds of nodal networks. We apply these statistics in particular to nodal 
patterns of random waves and of eigenfunctions of chaotic billiards. Of special interest 
is the observation that for monochromatic random waves, the number variance of the 
intersections with long straight segments grows like LlnL, as opposed to the linear 
growth predicted by the percolation model, which was successfully used to predict 
other long range nodal properties of that field. 



PACS numbers: 05.45.Mt, 05.40.-a, 03.65. Sq 
1. Introduction 

Entangled planar networks of convoluted lines appear quite often in various studies 
in physics and mathematics. To cite a few examples, recall the snapshots of polymer 
solutions, the level sets of rugged terrain, the trajectories of Brownian particles [Ij, 
the domain boundaries in magnetic materials or in simulations of random percolation 
[2] etc. To characterize such networks in a concise and quantitative way, one usually 
recourses to statistical measures, chosen to describe the specific properties of the network 
which are relevant to the problem at hand. Some studies give bounds on the length 
or the curvature of the nodal lines of wave functions in bounded domains, in other 
instances, the distribution of fractal dimensions of Brownian trajectories which provide 
an impression of their ruggedness [3j. Other examples are the distribution of the areas 
of connected domains in critical percolation which show a universal power law, or the 
(properly normalized) number of bulk nodal domains in billiard wave functions which 
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follow universal patterns that distinguish between chaotic and integrable billiards [4], 
[5]. None of these measures provide a complete description of the complexity of the 
network under study, and there is always room for introducing new measures which 
shed light on features which were not brought to the front by previous studies. 

In the present work we are interested in the point process generated by the 
intersections of a complex network of lines with a given reference curve. To the best 
of our knowledge, previous studies of the statistics of such point processes were not 
as detailed as the corresponding analysis for other systems (such as energy spectra, 
eigenvalues of random matrices, or one dimensional gas). In particular, while the density 
and correlation function were calculated, other physically meaningful observables which 
can be extracted from them, such as the number variance, rigidity and power spectrum 
were not investigated. 

Longuet-Higgins, in his analysis of the moving surface of the sea j6], derived the 
density and correlation functions of the zeros on a straight line. Berry and Dennis [7] 
studied the distribution of phase singularities in a random complex field. For a planar 
field, these singularities are actually the intersections of the nodal lines of the real part 
of the field with those of the imaginary part. In these cases (and others as well, including 
the present paper), restricting the field to a straight line allows using the formula derived 
by Rice [8j for calculating the exact correlation of the zeros from the known two point 
correlation of the amplitudes. Blum et al [4J considered billiard wave functions, and 
studied the number of intersections of their nodal lines with the billiard boundary. It 
was shown that the distributions of the (properly normalized) number of intersections 
distinguish between chaotic and integrable billiards. Johansson [Ij considered the 
trajectories of Brownian particles on a line, which start equally separated at t = and 
return to their original position at t = 2T without intersecting each other's path. The 
distribution of spacings at t = T, measured in units of the average spacing, was shown 
to be identical to the distribution of eigenvalues in the Gaussian Unitary Ensembles of 
N X N matrices (GUE). Baik and Rains [9] showed that for a certain class of discrete, 
non-intersecting random walkers, when not restricted to return to the original position, 
the limiting distribution of spacings matches that of the eigenvalues of the Gaussian 
Orthogonal ensemble (GOE). Such problems, when considered in (x,t) space, represent 
intersections of random planar paths with a reference line orthogonal to the t axis. 

In the present paper we shall study the points generated by the intersection of a 
reference curve with the nodal lines of three fields. 

(i) Random monochromatic waves [lOj in which are solutions of the wave equation 



Since no boundary conditions are required, the solutions are plane waves which 
propagate in arbitrary directions with wave vectors |fe| = k. The random wave 
ensemble can be constructed as a linear superpositions of plane waves: 



(1) 




(2) 
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(with r G M^, |fc^| = k and ^ 1), where kn/\k\ and 0^ are distributed 
uniformly and independently on the unit circle. Equivalently, one can also use 
random superpositions of solutions of the wave equation ([T]) in polar coordinates 

aoMkr) + 2 ^ aiJi{kr) cos{W + 0/) (3) 

with real coefficients a^, which are identically and independently distributed 
Gaussian variables, and where the phases are independent and uniformly 
distributed on [0, 27r]. 

The correlation function for this Gaussian ensemble is 

GMr,r') = .Mk\r-r'\) (4) 

(ii) Given a reference curve, we consider the normal derivative n • \/ip{r) of the random 
wave field ([2]), where n is the direction normal to the reference curve at r. For 
example, when the curve is the x axis, the two point correlation of this function is 

Gnrw(:^, x') = 2Ji{k\x - x'\)/{k\x - x'\) (5) 

The zeros of the normal derivative are not, strictly speaking, nodal intersections, but 
they do define a point process, which can be studied the same way. Furthermore, 
when the reference curve is straight, these zeros are the nodal intersections of an 
actual random field. 

The boundary modified random waves [H], were introduced in order to study the 
effect of a Dirichlet boundary on the statistics of the nodal set. They are defined 
in the upper half plane and are solutions of ([T]) subject to the boundary condition 

^^(x,?/ = 0) = (6) 

As explained in section [H the nodal intersections of that field with its boundary, the 
X axis, are identical with the zeros of the normal derivative of a Gaussian random 
wave field on this line. This fact makes the normally derived random wave field 
on a curved line a natural candidate for comparison with generalized boundary 
modified models, such as Wheeler's generalization for a circular boundary [12], or 
the semi-classical treatment of a general boundary by Urbina and Richter p^ . 

(iii) Random Gaussian fields with the short range correlation function 

GsRF(r-,r') = exp(-A;>-rf/4) (7) 

were introduced in to investigate the importance of the fact that the random 
waves correlation function decays very slowly. The correlation function above 
coincides with Grw at short ranges, but decays quickly at large values. 

According to the Uhlenbeck theorem [15j, the crossing probability of nodal lines 
of random waves vanishes. In [16j, the "avoidance" of nodal lines was defined, as a 
measure of the distance between nodal lines which avoid intersection in the vicinity of 
saddle points. The probability distribution P{d) of the avoidances d was computed and 
it displayed linear repulsion at small values of rf, P{d) oc d. However the coefficient of 
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proportionality was not the same as one would get from the application of Johansson's 
result to the case of Brownian trajectories. 

Having in mind the results of Johansson's paper on the one hand, and knowing 
that nodal lines avoidances repel linearly on the other hand, one might expect that the 
statistics of nodal intersections would also follow (at least to some extent) the predictions 
of random matrix theory. We shall show below that this is not the case. 

This paper is constructed as follows. In section O we discuss the derivation of the 
nodal intersection statistics from the known distribution of the amplitudes. First, we 
extract the necessary one-dimensional correlation functions which follow from restricting 
the field to a general reference curve (similar to the restriction of the amplitude to a 
straight line, which was done in [6j). Then, the statistics of the zeros is extracted from 
these one-dimensional amplitude correlations. An exact formula for computing zeros 
correlations from one dimensional amplitude correlations was derived by Rice in his 
analysis of random noise [17j,^. The derivation described here, uses a different method, 
and is similar to the derivation used by other authors ^18j,^. However, as opposed to 
these studies, we consider the case of a general curve, and derive the formulae for the 
line and the circle as special cases. Restriction to a general curve does not preserve 
the field's translational invariance. This makes the results more complicated, but the 
methods remain the same. Such a generalization is necessary in cases where the test 
curve is inherent in the problem we wish to model (such as the boundary of a quantum 
billiard), and not independent of the field. In [Appendix B| we describe the generalized 
calculation for the normally derived random waves. 

In section [3l we apply the above to the monochromatic random wave (Ej), using 
a circle and a straight line for reference curves. A universal correlation function is 
extracted for the semiclassical {k ^ 1) limit, and found to decay slower than the 
RMT correlations. The number variance and power spectrum are calculated and their 
special features, which are clearly different from the RMT predictions are shown as well. 
Over large ranges, the number variance is shown to grow as LlnL. This is in contrast 
with the linear growth we might expect from Bogomolny's percolation model [5], whose 
predictions for other asymptotic properties of the field (nodal count [4j and SLE driving 
force [19]) give a satisfactory match to numerical data. In section [H we discuss the 
zeros distribution of the normal derivative of a Gaussian field on the reference curve. 
In the framework described above this may be considered a different form of restricting 
the field to the curve. The semiclassical limit, number variance and form factor are 
calculated for this case too. 

Finally, in section [5] we compare the above results to those of other fields. The 
statistics mentioned above are calculated for the short-range Gaussian model ([7j) (which 
displays a different behaviour), and also evaluated numerically for chaotic billiards, 
which show a reasonable match to the monochromatic random wave model. The 
theoretical properties of the distributions considered are summarized in table [H and 
compared to the well known results for the Poisson process and RMT ensembles. 
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2. Statistics of zeros from amplitude correlations 

2.1. Nearest neighbour distribution 

When considering the statistics of a sequence of points, the nearest neighbour spacing 
distribution (where the spacing is measured in units of the mean spacing) is perhaps 
the most natural, and easy to evaluate experimentally. In figured] we show the nearest 
neighbour distributions for two point processes. The first is the intersections of the 
nodal lines of the random wave (RW) field ([2D with the reference curve. The second 
process (NRW) is the zeros of the normal derivative of the random waves on the curve. 
The distributions were generated by numerical simulations, and compared with the 
corresponding statistics of the random (Poisson) ensemble and the Wigner surmise 
(corresponding to the GOE ensemble) [20] , 



1.6 I 1 r 




s 



Figure 1: Density of the nearest neighbour level spacing. 

At short ranges, both ensembles display linear level repulsion similar to the GOE 
distribution, but with diflFerent slopes. The most conspicuous diflFerences appear at large 
spacings, and to get a clearer impression, we show in figure O the same data in semi-log 
plot. We observe two important differences. 

(i) On average, both distributions decay exponentially with approximately the same 
rate p{s) ^ exp(— 1.45). This decay is faster than the exp(— 5) Poisson decay, but 
slower the "Semi Poisson" [21j distribution, which decays like exp(— 2^). 

(ii) The overall exponential decays are decorated by persistent oscillations, which are 
clearly seen in the inset of figured The oscillations have slowly decaying amplitudes 
and their frequencies are ^ I.Att for RW and ^ 27r for NRW. 
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Figure 2: Nearest neighbour level spacing (log density). Inset shows persistent 
oscillations relative to the mean decaying curve. 

The results above demonstrate that the nodal intersection statistics are significantly 
different than those of other, well known, point processes. The purpose of this work is 
to investigate the reasons for the occurrence of these diflFerences. However, the nearest 
neighbour statistic is not easily amenable to analytic derivation. Rather, we use the two 
point correlation function, which is readily accessible, as will be shown in this paper. 
For the asymptotic 5 <C 1 limit, this statistic coincides with the nearest neighbour 
density. Furthermore, other statistics with well defined physical meaning, such as the 
power spectrum and number variance can be expressed in terms of the correlation. We 
will now give some definitions, and follow with an overview of a derivation for the two 
point correlation of nodal intersections with a general reference curve. 

2.2. Two point correlations 

A random field ip{r) is Gaussian if for any vector u whose elements 2/2, • • • , '^n are 
field amplitudes Ui = ip{ri) or partial derivatives of the field at some points ri, . . . , r^, 
the probability distribution of u is multivariate normal. The probability density of a 
multivariate normal random vector with mean value Uq is given by 



(8) 



where C is the covariance matrix Cij = Cov(t^^, Uj). For the random wave ensemble (Ej), 
Gaussianity follows from the n dimensional central limit theorem. 

Assuming smoothness of the field, means and covariances of the field's derivatives 
can be calculated by deriving the appropriate amplitude statistics. Combining that 
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with the fact that aU the elements of Uq and C are single point means and two 
point covariances, it follows that all statistical properties of a Gaussian field ip{r) are 
completely determined by its mean value {ip{r)) (which will be zero for the relevant 
cases investigated) and the two point covariance function G{r^r') = CoY{ip{r)^ip{r^)). 
As a consequence, we should be able to extract the distribution of the zeros from this 
function, and, in particular, the correlation function of the nodal intersections can be 
derived from the correlation function of the 1 dimensional restriction f{t) = ip{r{t)) of 
the investigated field to the given reference curve F : [0, L] ^ = r(t) {t is the natural 
curve parameter). 

The density and correlations of nodal intersections (NI) can be calculated directly 
once f{t) and its properties are known, in particular we shall use the two point 
correlations of the restricted field / and its derivative along the curve f{t) = d//dt: 

Co{t,t')^Cov{f{t)J{t')) 

C,{t,t') = Cov(/(t),/(t')) , Cr{t,t') = Cov(/(t),/(t')) (9) 

C2{t,t')^Coy[f{t),f{t')) 

We will also use lowercase c to denote the single point limit of these, i.e. Co(t) = Co(t, t) = 
Var(/(t)) etc. 

2.3. Restricting the field to a curve 

Before we write down expressions for Co, Ci and (^2, let us take note of special properties 
of these functions, which follow from the symmetry of our field. The fields considered 
are translationally invariant. This implies that the mean and the variance of the field 
are constant over the domain. It allows us to apply a constant linear transformation 
= Alp + 5, normalizing the field so that {ip{r)) = and Cq = {ip'^ (r)) = 1 for all r 
(this is true for equations ([2j) and ([3]), however we will keep Cq in our expressions below, 
to make explicit the effect of the multiplicative scaling). 

Furthermore, assuming smoothness of ip^ it follows that 

2(7/;(r)vV^(r)) = (r)) = 0. (10) 

In general, the restricted fields are not translationally invariant. Therefore the 
computation of their correlations requires special attention, and will involve the 
geometric properties of the reference curve, such as its curvature. However, in the limit 
t f , one can still project ffTOl) onto the curve, and conclude that Ci{t^ t') —> — Ci(t, f ), 
and in particular Ci = 0. 

The fields we consider here are also isotropic^, so the correlation function has the 
form ('0(r)'0(r^)) = G{\r — r'\). To make use of this property, we change the coordinates 

u = r + r' , d = r — r' , 

X An example of a random field which is not isotropic, is the non- monochromatic random wave models 
used in [6j, which have a finite average wavenumber vector (carrier wave) 
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so that derivatives with respect to u vanish. Derivatives with respect to d can be 
expressed in terms of the single parameter function G{d) (where d = \d\). 
This results in the following expressions for the correlation functions ([9]) 



Co = G{d) 

Ci = -r' ■ dC'id), Ci = r-dC\d) (11) 
C2 = -Y.^^ 



G" - ^ + ^5,,. 



did? d 



where r' = r(f ), and d = d/d. 



2.4' The nodal intersections statistics 
The density of zeros of / takes the form 

p(t) = Y,Kt-U)=5{f{t))\f{t)l where /(t,) = 0. (12) 

i 

FoUowing M. Kac [22], we use the Fourier representation: 

i //I ^ 

This integral can be solved using elementary properties of multinormal variables ([8]). 
For any multinormal variable x with zero mean, and real vector fc, the following identity 
holds: 

^-^ ^ kiCi^jk^ where C^j = Cov(xi, Xj) (13) 

Using the Gaussianity of the field /, this identity can be applied to the mean density. 
The resulting integral separates into simple Gaussian integrals, and we get the following 
result for the density: 

(P) = -\{f (14) 

We now move on to two point statistics. When dealing with "point processes" (such 
as the zeros of /), whose amplitude is a sum of delta functions, {pp') normally contains 
a delta function dit t = t\ corresponding to the strong autocorrelation. It is customary 
to subtract this term, to regularize the correlation function at t = f. Thus, the two 
point correlation function is defined (e.g. in [20]) as 



(e'^'^) = exp 



= {pmt^))-s{t-t^){p{t)). 

Using the same approach as above, we write: 
W> = (i^ / ^^^^f^-'^'d - e'")e'"'(l - e-'/')) (16) 
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To make use of (1131) . we write down the covariance matrix of the four variables 
(/) /') /) /') appearing in this expression (see definitions in (jHl)): 

/ Co Co Ci\ 

(a r\ 

(16) 



M 



Co Co Ci 

Ci C2 C2 



\Ci C2 C2 y 




{A,B,C stand for the 2x2 submatrices) . 

Applying f|T3l) to f|T5l) we get more complicated integrals, which are nevertheless 



solvable (derivation in [Appendix A[). The resulting expression for R{t,t') is 



R{t,t') = — 1^13/2 (^Vl - + carcsin(c) 



^2 |^|3/2 

where a = C2\A\ - cq\C\, c = {C2\A\ - Cq\C\) /a, \A\ = cq^ - Cq^ and \C\ = -CiCi (here 
1^41 and \C\ stand for determinants — not absolute value). 

From this, the normalized correlation coefficient TZ = R/ (p)'^ — 1^ from which we 
will derive the other two point statistics, follows immediately: 

n{t,t') = (yr^ + carcsin(c)) - 1 (17) 

3. Simple curves in monochromatic random fields 

The monochromatic Gaussian field ([2j) is often used as a statistical model for 
eigenfunctions of chaotic billiards. In this context, the "semiclassical" regime (large 
k) has special significance. In terms of our statistical system, this means that we would 
be particularly interested in the k ^ n limit (where is the curvature of the reference 
curve). In this regime, the first approximation for the curve is a straight line, and the 
second one is a circular arc. These two cases are also important for their simplicity — in 
both cases, the domain distance d between two points on the curve depends only on the 
curve distance |t — f |, so ffT2l) assumes a simpler form. 

For the straight line, r = d is constant and = \t — t'\. We get: 

C^ = G{d), Ci = -G\d), C2 = -G"{d) (18a) 

For a circle of radius r, define a = {t — f)/(2r), so d = 2rsina, and by simple 
geometrical identities, we get: 

C^ = G{d), Ci = - cos aG\d), 

C2 = sin' - cos' aG'{d) (186) 
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3.1. Normalized correlation for the line and the circle 

Substituting (j4j) in (I18aj) and fll8&j) . we get the amplitude correlations. The results for 
the circle are (the straight line is given by their a ^ limit): 



Co = Mkd) 

Ci = kcos{a)Ji{kd) (19) 
C2 = k^(Mkd)cos\a) ^'^^"^^ 



kd 

For the density, we substitute C2 = ^2(0) = in ffT4l) getting (p) = k/{^/27r). We 
use this to select a curve parameter which measures length in units of average spacing 
(so we can expect meaningful results at the (p) 00 limit): 

s = {t' - t){p) = {t- t')k/{V2n). 

In the new units, a = s^/2ti / {2kr) and the argument of the Bessel functions is 

— - — > V2tis 

2kr I kr^oo 

For the NI normalized correlation, we substitute ffTOl) in ffTTI) . getting: 

[Jo{9) cos^ a - Ji{9)/9] [l - Jo{9f] - Jo{9)Ji\9) cos^ a 
^~ [l-Jo\9)] /2-Ji\9) cos^a 
^ 1 - Jo^(^) - 2 Ji^(e) cos^ a f /- . . A , 

7^= ^^^7^ Vl-c2 + carcsinc - 1 20 

[l-J.\9)f" ^ ^ 

To get the universal distribution [k independent, high density limit), we take 
k ^ 00 keeping s constant. 

[Me) - MO)/e] [1 - Mof] - ue)j,\e) 



c = 



[i-Jo\e)] /2-j,\e) 



1 - j,\e) - 2j^{e) ( ^ — - . ^ 

iZ = — V 1 — + c arcsm c 



Note that in the case of the straight line (/^ = 1/r ^ 0), we have 9 = \/2i\s and a = 
exactly, so there is no k dependence, and the universal function fl2Tl) is exact for any 
k (this is expected, because with no boundary and no curvature, the system has no 
natural scale). 

To further study the universal correlation function fl2T]) . we derive its asymptotic 
expansion. The first few terms in the expansion of V^i^s) near 5 = are: 

This rises linearly from —1, with a smaller slope than the GOE ensemble (other RMT 
ensembles are not linear at this limit), however, as is shown in figure [3l it rises quicker 
to positive values, resulting in a smaller range of high repulsion. 
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Figure 3: Normalized Correlation: NI of RW vs RMT level spacings. 



For large values of 5, we find decaying terms, oscillating in frequencies which are 
multiples of the dominant frequency cu = 27^^/2. The asymptotic expansion has the 
form: 

3 ^ n 

n{s) ^ Re V V <7.,^(i7r)e^-(— + 0(5"^) (23) 

where qn,m{,x) are the foUowing polynomials: 

qifi = 1, = 9 

_ 25 _ 49 39 _ 121 

^2,0 — Q2,l ~ "3" + ~2^' ^^2,2 — -J^ 

169 17 2 1369 511 205 2 

_ 1681 55 _ 529 

g3,2--^ + Y^x, q3,3-^ 

Writing the leading terms explicitly, we have 

7^ ^(l + 9sin(6c;5)) (24) 

TTUJS 

This decay is slower than the decay of random matrix ensembles for GOE 
and GUE, and x [oscillating part] for GSE). The slow decay has a considerable effect 
on all statistics that probe large distances, as will be shown in section [321 

When we look at intersections on a full circle, s goes over all values, from —kr/\/2 
to kr/^/2. This range contains values large enough to make fl2TI) unsuitable for 
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approximating the correlations. The leading terms in the asymptotic expansion of fl20l) 
(the equivalent of fl24l) ) are found to be: 

n{s) = {As + Bs sin(a5)) + 0(^-2) (25) 

where 



sinc| 




(with normalized sinc(x) = sin(7rx)/(7rx)), are slowly varying functions, starting from 
Aq = I^Bq = 9 and Qq = uj (as expected from fl24l) . which holds for s <C fcr), and ending 
in A±s = ^^B±s = 1 and f^^s- = 4:^/2 for aS = kr/^/2. The functions 74^, 5^ and Qs 
are illustrated in figure [H As we increase fc, the rate of change becomes slower, but the 
normalized length of the curve increases, and the value of the coefficients for maximally 
separated points (1^1 = kr/^/2) remains the same. 



/ 2iis 

2 + cos 




2 - ^ . . 

°0.0 0.2 0.4 0.6 0.8 1.0 

X 



Figure 4: Asymptotic coeflBcients and frequency for RW NI on half a circle as a function 
of X = ^/2s/ (kr). 



Before proceeding to compute other statistics, we would like to point out that the 
oscillations in p{s) (see section 12.11) have a frequency which is numerically consistent 
with half the dominant frequency uu = 2i{\f2. 



3.2. Number variance and power spectrum 

The variance of the number of intersections in a finite segment of length L is often used 
as a measure for the "rigidity" of the distribution (a completely rigid set of points will 
have zero variance). The power spectrum [17j measures periodicity in the distribution 
(which is also manifested by oscillations in the correlation function). These statistics 
are given by integral transforms of TZ which are hard to solve directly. Instead, the 
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asymptotic expansions of the correlation function can be used to expand them in series 



at regions of interest, as explained in Appendix C 



The variance for the number of intersections on a finite segment of length L can be 
calculated from the normalized correlation by [20j : 

S2(L) = L + 2 [ (L - s)n{s)ds 



(26) 

^0 

We can use f l22l) and fl23l) to find the asymptotic behaviour of this statistics. For small 
L, we find 

S2 ^ L - + -L^ + ^^L' + 

48 23040 19440^2 

The first two terms follow from the formal definition of S^, and therefore coincide with 
the corresponding terms in the RMT ensembles. The third term corresponds to the 
probability of finding two intersections in the small interval [0, L]. It is of order L^, as 
in GOE (again with a smaller factor). 
For large L, we find 

- ^(L In L - L) + L(l + 2M^w) 



+ 



-^(lnL + l)-2Mj, 
1 



RW 



^ Re[g3,o(i7r) + 2i7rVe'"T +0(L-2) 

[nuj) L 

where and A^kw are fixed constants, given by 

1 



(27) 



n{s)ds + 



n{s) 



TTUS 



ds 



-0.336 



TTUJ 



COSCJ 



1 



s7l{s)ds 







+ 



1 25 

sTZ(s) (1 + Qsincjs) ^ 

. ^ TTCJ^ ^ 4(7rcj)'5 

0.0826 



ds 



Writing down the first terms of fl27l) explicitly, we have: 



E2(L) 



~ — L In L + 

TTiO 



+ 2M 



TTLJ 





RW 



25 



(27ra;)^ 



lnL-2 ( A^^w + 



L 
25 



(27ra;)' 



(28) 



Note that this variance is asymptotically larger than that of the random (Poisson) 
distribution (where H'^ ^ L). However, due to the numerical constants, the LlnL term 
dominates only when L > exp(7r6j(A^^^ + 1/2) — 1) 36, and only at L 32 x 10^ 
does it become larger than the Poisson variance. At segments whose length is a few 
average spacings, the variance already becomes significantly larger than the "rigid" 
RMT ensembles. In figure O the variances of several distributions are compared. 
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Figure 5: Number variance in segments of length L. For large L, the NRW curve is 
linear, hence parallel to the Poisson curve. The RW curve rises faster (the crossing point 
is at L 32 X 10^, beyond the range of this plot). 



It should be noted that in the case of independent site percolation (for grids with 
bounded, convex cells), a straight line passes only once through each cell. Hence the 
sign changes along the line form a Poisson process, with linear number variance. The 
fact that the asymptotic behaviour of the random wave model behaves differently is of 
interest, because the predictions of the percolation model [5j matched the numerical 
data for other parameters which characterize the behaviour of RW nodal lines in the 
high energy limit. Namely, the nodal count (number of nodal domains contained inside 
a given area) [4j and the SLE driving force [19j. 

The form factor (scaled power spectrum) for the intersections is the Fourier 
transform of the scaled correlation function (up to the subtraction of a 6{r) term). 



Kir) 



^2^Ts {P{^)P{X + S)) 



Sir) 



(29) 



In figure [6] this is plotted and compared to the corresponding RMT values. 

By f l24l) it is clear that K diverges for r = (since I/^ttujs) is non integrable). 
Similarly, from the form of equation fl23l) , we expect singularities at integer multiples of 
^/2 (corresponding to the angular frequency uj = 27r\/2). To quantify this, we use the 
asymptotic expansion of TZ to expand K{t) in small regions around the singular points 
T = V2n + 5 for 5 < 1. 
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For the divergence at (r <C 1), we get 
2[7 + ln(2^)] 2 



K{r) - 1 + 2A^^ 



RW 



Inr 



TTLU TIUJ 

This is a logarithmic divergence. It drops down to 1 at tq ^ 8 x 10~^ (which is too small 
to be observable in figure [6]). The parameter K{0) = 1 + J^7l{s)ds can be viewed as a 
measure of asymptotic spectral rigidity [23j. In our case, we see that we have "infinite 
softness" — the slow decay of the correlations outweighs the short distance repulsion. 

One notable feature of figure O is the finite jump at r = ^/2. Expanding for this 
region, we find 



K(y2 + 5^ ^ A + B5 



39 ,„ 49 



where A 



1.25 and S ^ 1.53 are numerical constants. 
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Figure 6: Form factor for RW Nodal Intersections. The logarithmic divergence to +oc 
at r = is too steep to be observable in this scale. 

For the case where the reference line is a circle (inset of figure [6]) , we observe strong 
oscillations in the form factor for r < \/2. These oscillations are caused by the "drifting" 
of the frequency in (ESI). Their frequency increases for large values of fc, and we would 
expect them to cancel out when averaged over a wide enough window of frequencies. 



4. The normal derivative of a Gaussian field 



When dealing with solutions of boundary value problems (e.g. two dimensional quantum 
billiards), the boundary of the domain is a curve that plays a special role. It is therefore 
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worthwhile to investigate its intersections with the nodal lines. For a Dirichlet boundary 
condition, the boundary is part of the nodal set, and the derivative of the eigenfunction 
in the tangent direction is zero. At the points where another nodal line intersects the 
boundary, the derivative in the direction of the other line will also vanish. Since nodal 
lines of solutions of the Helmholtz equation intersect at right angles, the intersections 
we wish to explore are also zeros of the normal derivative of the field with respect to 
the curve. This fact makes the normal derivative of other fields, such as the unbound 
random wave, natural candidates for comparison with the statistics of such intersections. 

If fact, the normal derivative of a Gaussian random wave field is directly related 
to the "boundary modified" field mentioned in section [H The boundary modified 
random waves [llj, with correlation function approaching (ED, might be realized by 
symmetrization of the unbound random wave field ([2j): '0brw(^5 v) = ^0(^5 v) —y) 
(where ip = t/^rw of (IID). In this case, the correlations on the boundary itself can be 
defined as the limit of correlations on a line parallel to the x axis at y ^ 0+. To get 
a meaningful limit we must introduce a ^-dependent scaling, to keep (V^brw^) constant 
as we approach y = 0. Denoting (l)y{x) = '^brwI^,?/), the average square amplitude for 
a line at height y is: 

{<i)y{xf) = y) - ^(x, -y)f) = 2 {Co- Co{2y)) 

The correlation function Co is symmetric, so in the generic case (and, in our case 
specifically), its value for small d can be approximated by a square function C{){d) ^ 
Co + \c"()d^ . Thus, {(j)y{xY) ^ (2^)^|c^^o| (note that c^'o is negative), and the scaled field 
will be: 

(py{x) = ^ 

^ 1 "ipjx.y) - ^Jj{x, -y) ^ 1 d%lj{x,y) 

^ ^A^\ ^A^\ 9y 

Therefore, the properly scaled boundary modified field is proportional to the normal 
derivative of the random wave field used to generate it. 

Jf..!. Amplitude correlations of the normal derivative on the reference curve 

The normal derivative of the field across the curve is given by g{t) = n{t) • V'0(r(t)), 
where n{t) is the unit normal vector of the curve at point t. Repeating the computation 
described in section 12.31 we can calculate its correlation functions. 

{g{t)g{t)) =^n. n, 



G"{d) 



G'{d)\ didj ^ G'jd) ^ 



d 



d^ 



d 



n 



(30) 
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For the circle (as before a = {t — t^)/(2r), d = 2rsinQ{), we have the identities 
n ' n' = cos(2a) and • d = —n • d = sin a. Inserting these in fl30l) we get: 



Co = {gg') = - cos' + sin' aG%d) (31a) 

til 

The other two correlation functions are calculated in a similar manner (for a general 
curve) in [Appendix B[ For the case where the curve is a circle, they assume the following 



forn|§: 



Ci {gg') = cos a (l — 3 sin' a) f — j 



— sin' a cos a 

|^G(3)+4^) (316) 
C2 = {gg') = - 2 (cos(2a) + ^ sm2(2a)) 1 {g" - ^ 
l--sm^(2«)j^ 

- sin^ a (^os^ aC^^^ + 4 cos(2tt) ^) (31c) 

The correlations for the straight hne are easily derived from equations fl31all - fl31c|) 
by taking the limit a ^ 0: 



(32a) 



G' 




~~d 




G" 


G' 


T ~ 






2 


d 





= V - ^ ( - ? ) (32c) 



^1 = - - :^ (32&) 

d 

To study the normally derived field corresponding to the monochromatic random 
ensemble of section [3l we apply the results above to the specific correlation function, and 
discuss the resulting statistics. For simplicity, we discuss the case of a straight line first, 
and follow with comments about the main diflFerences in the case of a circular curve. 

For the straight line, the amplitude correlations are derived from fl32al) - fl32d) . 
To get a standard scaling with cq = 1, we first multiply the field by a/2/A:, getting 
G = 2Jo{kd)/k'^^ and the following correlations for the normal derivative: 

^ ^ 2Ji(fcrf) 
^ kd 

§ This form was chosen to make the a ^ hmit more evident. The forms in [Appendix B are more 
suitable for calculations involving Bessel functions or generic derivatives. 
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2J2{kd) 



kd 



2 2 f J2{kd) 



Ci = k- 

4^2. Nodal intersection statistics for a straight line 

We proceed as in section [3l First, the density of intersection is calculated using ffT4l) : 
1 [c^ k 

We use this to rescale the distance to units of average spacing s = [t — t')k / (27r), so the 
argument of the Bessel functions in fl33l) becomes 9 = 2tis. Substituting these results 
in (ITTI) . we get: 



c = 



9^ 9^ - W - 16 



2 



^ = _ ^^^2)3/2 ( V 1 - c2 - c arcsm c j - 1 (34) 

The asymptotic expansion of TZ for small s is: 

+ + + 324^ 

As before, this grows linearly from —1. The slope is larger than that of the unmodified 
RW, but smaller than the slope of the GOE correlations. As shown in figure [71 it rises 
quickly to positive values. 

The expansion for large s is: 

2 —555 
^(^) ^ ^^-^ - 25sin(4^5)) + 4 cos(4^5) 

1 / 69 12789 . , . 

+ sm(47r5) . (35) 



27r(27r5)^ V 2 

7^.(5) decays as which is faster than the decay of the correlations for the 
corresponding RMT ensembles. As in fl23l) . each asymptotic term has oscillating parts, 
with frequencies which are multiples of the dominant frequency cjnrw = 47r. As 
in section 13.11 we find that the oscillations in the nearest neighbour density of the 
NRW, described in section 12.11 have a frequency which is numerically consistent with 
half the dominant frequency uu = An. 

With the decay of the correlations, both J TZ{s)ds and J sTZ{s)ds are finite, 
so from fl26l) it is evident that should increase linearly for large values of L. The 
asymptotic expansion calculated from fl35l) is: 

^ (1 + 2M° Rw)^ - 2Mkw (36) 

with 



POO 

•Mnrw= / 7^(s)ds ~ -0.2582 
Jo 
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Figure 7: Normalized Correlation: NI of normally derived RW. 



and 



POO 

-^NRW = / sn{s)ds 0.00617. 

Jo 



As seen in figure El the asymptotic variance is larger than RMT ensembles, but smaller 
than Poisson. 

Much in the same way, from fl29l) . it is evident that the form factor will now be 
finite at r = 0, and will be continuously different iable for r > 0. Irregularities appear 
in second and higher derivatives at multiples of the dominant frequency cjnrw (i-e. 
Tn = na;NRw/(27r) = 2n). The asymptotic expansion for r <C 1 is given by: 

_9 

—2 



K{r)^{l + 2M'^^^) + -THnr 



A(7-^ + ln27r)-(27rfA^^RW 



with 



>^NRW = ^ s^n{s)ds + ^ (^^(^) - ^ 0.00579 

We now have a finite K{0) ^ 0.4836. This value is smaller than the K = 1 oi the 
Poisson distribution, but larger (softer) than the RMT ensembles, for which K{0) = 0. 



4-3. Statistics on a circular arc 

The case where the reference curve is a circle can be treated similarly. The amplitude 
correlations of the monochromatic waves (corrected form for fl33l) ) can be calculated 
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Figure 8: Form factor for the NI of normally derived RW. 



directly from (IB. 51) . It should be noted that the amplitude variance Cq = Co(0) remains 
^/c^ (as for the straight line). However, for any finite fc, the variance of the derivative 
(C2) is different than the value obtained for the straight line, and induces a larger density 
of intersections. 

The corrected density of intersections for the circle is: 



P = 



k_ 
2^ 



271 



k_ 

2^ 



1 + 2 



The 0(/^^) increase in density, which is due to curvature, affects the choice of the scaled 
curve parameter s = p - {t — f) and thus all the correlation formulae. It is convenient to 
express the resulting statistics in terms of the dimensionless parameter 

which goes to 1 when kr ^ 00. With this, the small s expansion of the normalized 
correlation is: 



7^ 



1 + 



+ — /?'(-24 + 357/3' - 820/?^ + 500/3^)5^ 
576 

For large values of s (the circle is finite, and \s\ can only get values up to 
l^r = (1 — Z?^)""^^^), the constant dominant frequency lo = An of fl36l) is replaced with 
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a drifting frequency (changing very slowly with s) 
f^s = 47r/3sinc^^^ , 



similar to the situation in section ISTTl However, in the case of the normally derived field, 
there is another notable difference between the circular and the straight curve. Namely, 
for any finite value of fcr, we get 0(5~^) and 0{s~'^) terms, which do not appear in the 
case of a straight line. The leading term of the normalized correlation for 5 ^ 1 is 

4 sin^ a[(4/?2 cos^ a-1? - (4/?^ cos^ a + 1)^ sinfa^)] • 

TrilsS 

where a = 7rs/{kr) is slowly varying up to 7r/2. When we are in the 1 <C 5 <C fcr region, 
a is very small, and this term is 0{a^). Hence, it vanishes in the straight line limit. 
The 0(5"^) term, contains the frequencies 0, and 2Qs- In the low a region, it is 

(again, vanishes for the straight line). The 0(5 ^) term is the 
first one that does not vanish when a ^ 0. 

The effect of these terms is manifested by the number variance statistic. In figure O 
the number variance is plotted for three values of kr. As long as L <C fcr, the variance 
stays close to the values corresponding to a straight curve given by fl36l) . However, 



when the angle covered by the segment becomes non negligible, 0{a^s~^) and 0(a^5~^) 
terms of the correlation become important, and the variance increases above the linear 
asymptote. 



400 
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Figure 9: Number variance for zeros of NRW on a circular curve, for three different 
values of kr. The length of the segment is bound by the total length of the circle, 
marked by circular dots on the plot. 
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5. Comparison with other fields 

The statistics considered above can be evaluated numerically for the nodal intersections 
of any scalar field with a given reference curve. If a random Gaussian model is given 
for the field, the techniques demonstrated above can be used to extract the theoretical 
functions predicted by the model, and those could be compared to the numerical results. 
This gives a useful tool for verifying statistical models for such systems. For the case 
of chaotic billiards, the numerical results can be compared to the predictions of random 
wave models. Also, we can use these statistics to compare other random wave models 
with the monochromatic random waves. In ^4j, the field with correlation ([71) was 
compared to the random waves model with respect to special statistics proposed there 
(which were related to the probability of finding a nodal line inside a reference tube). 
The nodal intersection statistics provide another way to compare these models. 



5.1. The Gaussian short range field 

The short range field ([71) was chosen so that its short distance behaviour will be the 
same as that of the monochromatic random wave field, so we expect the average level 
spacing to be v^tt/Zc. Substituting the correlation function ([71) in fll8&l) . and denoting 
G exp(-6>^/4), 9 kd, we get: 



Go = G, Gi = -kecos{a)G, 

G2 = h^[co8{2a) - ^e^co8\a)]G 

(so 02 = /c^/2 as expected, and 9 
for the straight line becomes: 

1 



\/27r5 as before). The nodal intersections correlation 



e ^ 



e 2^ 



1 



e •2" 



nis) 



1 — e 2^ 



3 



Vl — + carcsin(c) ) — 1 



(l-e-2^ , 
The expansion for small s is: 



(37) 



54 



and for large s, it is: 



iws)^ -2{7rsf + 1 



As can be seen in figure llOal the slope at is higher than that of the GOE, and the 
correlations decay to faster than all other distributions considered here. Thus, this 
field has smaller repulsion and it may be considered the closest to the random Poisson 
distribution. 



Nodal Intersections With a Reference Curve 



23 




(a) Correlations (b) Form Factor 



Figure 10: SRF nodal intersections compared to other fields. 



The number variance and form factor also refiect this property. Denoting 

POO 

^SRF = / T^{s)ds - -0.2141 

POO 

•^SRF = / sn{s)ds 0.02902 

Jo 

roo 

-MsRF = / s^n{s)ds - -0.00434 

We find that the asymptotic number variance is linear with coefficient (1 + 
^ 0.5718, which is larger than the normally derived random waves (thus closer 
to Poisson). 

The form factor has a smooth form (because fl37l) is not oscillatory). The expansion 
for small r is: 

KiT)^il + 2Ml^^)-A7r'Ml^^T' 

We get a relatively large K{0) ^ 0.5718, and a quick approach to the limiting value of 
1 (which is consistent with the notion that this is closer to the Poisson distribution). 



5.2. Chaotic billiards 

Eigenfunctions of chaotic billiards provide an ensemble of functions with complex 
behaviour. In the semiclassical limit, they are well modelled by Gaussian random 
wave ensembles (the monochromatic random wave ensemble is used to model the 
eigenfunctions in regions which are far enough from the boundary, and various boundary 
modified ensembles model the field close to the boundary). 

The statistics considered in this paper were evaluated numerically for eigenfunctions 
of desymmetrized Bunimovich stadium (0.5 x 1 rectangle joined with a quarter of a circle 
of radius 1) and desymmetrized Sinai billiard (1.2 x 1 rectangle, with a quarter of a circle 
of radius 0.5 cut out of one of the corners) — both with Dirichlet boundary conditions. 
For the Stadium, 1500 eigenfunctions were taken, with wave numbers ranging from 
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fc=110to/c = 165. For the Sinai billiard, we used 10000 wave functions, with wave 
numbers from 350 to 500. 




(a) Stadium billiard (b) Sinai billiard 



Figure 11: Internal reference curves in the quantum billiards 

Fixed reference curves were chosen in the interior of the billiards, as shown 
in figure [m For each wave function, the sequence of intersections of the nodal lines with 
these curves was calculated, and normalized to unit average spacing according to the 
corresponding wavenumber. As shown in figure [121 the nearest neighbour distribution 
agrees very well with the predictions of the random waves model. The form factor is 
also very close to the expected curve. 




(a) NN spacing (b) Form factor 



Figure 12: Billiards: NI with internal line. 

The number variance of nodal intersections on segments of the reference curves 
is plotted in figure [131 For short segments, the results are close to the predictions 
of the random waves model, but deviate considerably from this model as we move to 
segments of large normalized length, especially in the case of the Sinai billiard. This 
deviation happens when the (unsealed) length of the segment is comparable to the 
billiard dimensions. On such scales, it seems likely that the geometrical details of the 
billiard would have an effect on the statistical properties of the wave functions. 
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Figure 13: Number variance of NI with segments inside chaotic biUiards. 
6. Summary 

The analytic results obtained for the three Gaussian fields discussed in this paper are 
summarized in table [U and compared to their analogues in the well known random 
matrix and Poisson ensembles. For short distances, the normalized correlation function 
for all the distributions considered (except of the Poisson distribution) approaches —1, 
so the leading term of 1 + TZ{s) is listed, for 5 = e <C 1. Similarly, the number variance 
for short intervals is asymptotically L — (again, except of the Poisson case), so the 
leading term of R2 = S^(L) — L + = E[n{n — 1)] is listed, (for L = e <C 1, this is also 
twice the probability for finding two points in the interval). For large distances, we list 
the leading terms of and TZ. For the form factor, we list the "asymptotic rigidity" 
K{0) (which is directly related to and A in large intervals [20j,[23j), the singular 
frequencies where K{t) or its derivatives diverge and the type of divergence (number of 
the first diverging derivative) in the first non zero singular frequency. 
We have seen that 

• Nodal intersections (NI) with a reference curve can be used as a statistical means 
to obtain insight into the properties of complex networks of lines. The NI statistics 
of the random wave (RW) model display characteristic behaviour, which is quite 
different from the statistics of other point processes. Eigenfunctions of chaotic 
billiards match these statistics well, with some differences on scales comparable to 
the dimensions of the billiard. 

• For the RW and some other Gaussian models, the two point correlations and other 
statistics related to it can be derived by analytical means. 

• The nearest neighbour distribution is harder to obtain (since it depends on n-point 
correlations with n > 3). However, numerical evidence suggest that the frequency of 
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Correlation Form Factor Number Variance 



Process 





1 + 7^(e) 


n{s > 1) 


K{0) 


Divergences^ 


^2(e) 


S2(L > 1) 


roisson 


1 


u 


1 




2 


T 

L 


Oxlr 


Z.4D / 6 


4o. / 5 e ^ ^ 


u.o / z 




u.ozze 


u.o / Zi^ 


NRW 


1.234e 


0.023S-3 


0.484 


(3){2,4,...} 


0.411e3 


0.484L 


RW 


0.617e 


0.036S-1 


oo 




. .} 0.206e3 


0.072L logL 


GOE 


1.645e 







(4)|1} 


0.548e3 


0.203 log L 


GUE 


3.290e2 


-0.051S-2 







0.548e4 


0.101 log L 


GSE 


11.55e4 


0.25s~-^ cos{ujs) 





(0){1,2} 


0.770e*5 


0.051 log L 



Nonzero values of r, for which K{t) or one of its derivatives diverge. The number in 
parentheses specifies the diverging derivative at the first of these points. 



oscillations is related to the dominant frequencies calculated for the corresponding 
two point correlation functions. 
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Appendix A. Derivation of the zeros correlation from amplitude statistics 

In section 12.41 the correlation of the zeros of a Gaussian random function on a curve 
was expressed in terms of the Fourier integral ffTSl) : 

/ /\ / 1 / dC dry d^' d^ .^^ ■ f .^;f; ^ \ 
(PP) = ^^ J e^^(l-e^Oe^Ml-e^M^ 

We note that the integrand contains four terms of the form exp{iVi • /), with 

Vi e {(e, r, 0, 0), (e, r, V, 0), (e, 0, v'), v, V)}, 

and f = {fJ'JJ'). 
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If we use f lT3l) to evaluate the statistical average, these terms are replaced by 
exp(^—^Vi^ Mvij ^ where M is the covariance matrix of /, given in ffTGl) . To continue, we 
first integrate the ^ and coordinates. To separate them out, we use a simple block 
matrix identity. 

For any pair of column vectors x^y ^ and three n x n matrices A^B^C ^ where 
A and B are symmetric and A is non-singular, the following identity can be simply 
derived: 

V) 5 ) ( ^ ) = + ^")^^^ + ^) + y"^y ^^-1) 

where B = B- C^A'^C and y = A'^Cy. 

In addition, using the same notation, one can also deduce the determinant identity 

= detAdetB, (A.2) 

JD I 

which we shall use below. 

Applying flA.ll) to the four terms, with = (x^^yi^) (i.e. x^ = (<^,<^0 
yi^ G {(0, 0), (ry, 0), (0, 7y^), (ry, 7y^)}), the x dependent integration for each of the four 
terms reduces to the usual Gaussian form: 




/ 



VdetA 

This does not depend on y and therefore contributes a constant multiplicative factor. 
The remaining integral assumes the following form: 



oo 



(using a, 6 and c for the matrix elements of B). 
This can be rewritten as Ii + /2, where 

A=/r^(l-e-'-)(l-e-"'=)=2.vS 

J J — OO '/ '/ 

and 

h = [T ^e-'^^'^^'+'^''\e-^^^' - 1) 

J J — OO V V 

I2 is not well defined, but taking the Cauchy Principal Value for the integrals, anti- 
symmetric integrands do not contribute. After symmetrization and change of integration 
variables, we get: 



Vd) ir [cosh(cx2/) - 1] 

J j-00 ^ y 



with c = c/^/ab. One way to solve this integral is to derive the expression by c 
twice, and integrate back after the x and y integrations. Another way is to express 
f{xy) = {xy)~'^ [co8h{cxy) — 1] as a power series in (x^)^, and integrate term by term 
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(each term becomes a square of a simple Gaussian moment). Either way, for any |c| < 1, 
we get 

I2 = 27r^/ab —1 + Vl — c^ + c arcsin c , 

To show that the convergence condition |c| < 1 holds, we note that 

|c| < 1 <^ det B = ab - > 0. 

However, from flA.2[) we have det M = det A det B. Since both A and M are covariance 
matrices, they are positive definite, and their determinant must be positive for any 
non-degenerate case. This can only be consistent if det B is positive too. 

Inserting Ii and I2 into the expression for the normalized correlation we find ffTTI) 

TZitA') = — ^.^ f a/1 — + carcsin(c) ) — 1 

C2{deiAf'^ V V 

Appendix B. Calculating the amplitude correlations for the normal 
derivative of a Gaussian field 



When the investigated field is translationally invariant and isotropic, the correlation 
functions Co, Ci and C2 for the "normally derived" function g{t) = n{t) • V'0(r(t)) can 
be calculated in terms of the curve parameters and reduced correlation function G{d). 
The derivation is straightforward but tedious, and follows the same mechanism used 
in section [2l3l and section [4Tl In this section we recite the main steps and results. 

As in previous sections, we start by using linearity to move the derivatives out of 
the statistical average, getting fl30l) for the first correlation function Cq. This expression 
involves second order partial derivatives of the basic correlation {ipijj^)^ so the expressions 
for Ci and C2 involve third and fourth order partial derivatives correspondingly. 

Using shorthand notation where tagged parameters refer to values at f^di = d/dri 
and d'i = d/dr'i^ we have 

Co= Y.nin',did',G{\r-v'\) 
Ci= Y^d'trnn'.did'.G 

C2= Y,^t&tn^n',^0,G (B.l) 

The derivatives of n in these expressions can be evaluated using Frenet formula 
dtTi = —nr (where n is the curvature at t), and the partial derivatives of G can be 
expressed in terms of derivatives with respect to d {G' ^G" ^G^^\G^'^^) and the distance 
vector d. 

First, we calculate the partial derivatives of G using the following rules for 
derivatives of the distance d=\d\ and direction vector d = d/d (where d = r — r') . 
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the 



To express the results in a manifestly symmetric way, we define <^ii,^2v,^iv' 
dimensional symmetric direction tensor of degeneracy fc" as the symmetrized product 
of k Kronecker deltas and N — 2k direction vectors d. Specifically: 



didj 7 



d} 



Sijdk + Sikdj + 5jkdi 



Sijdkdi + Sikdj di + 5iidjdk + 
Sjkdidi + 5jididk + Skididj 

^fjkl ^ ^ij^kl + ^^z/c^^j7 + SiiSjk 

We also introduce the notation for the specific combinations of derivatives of 
G up to order n, which satisfy the following recurrence relation^: 



Specifically, we have: 



1 



- -G^ 
d 



G[4] = Q(A) _ 6^(3) + 



15 



G' 



With these notations, the partial derivatives of G are given by: 
diG = GWdf 



d 



d^dMG = + ^G[3i4i + 4^^145 



Substituting these in fIB.ll) . we get 

Co = J]-n,n',(GPl4? + ^G[^l4.?) 



'ijkl 



(B.2) 



which is equivalent to fl30l) . Similarly, 

Ci = ^ n'nir'jdidjG + ^ nin' jr'j^didjdkG 



(nd)(r^d)G'^' + {nr') 

C6 



II These relations are almost identical to the Bessel recurrence relations. Hence, for the case where 
G{d) = Mkd), we get GW = (-fc)" J„(fcd) 
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(B.3) 



C2 = - nn'rir'-dijC + dirjr'j^dijkG + nin' jrkr'idijkiG 



— 



{rd){r'd)G^^^ + {rr') 



+ [(dr)(r'<l) + {dr'){fd) + (dd)(rr') 

+ (nd)(n'd)(rd)(r'd)GW 

+ \(nd){n'd)(rr') + (nd) (r'd) (rn') 



+ {n'd) (rd) {nr') + (rd) (r'd) {nn') 
+ [inn'){rr') + {nr')irn')] ^Ct^l 



J d 



(B.4) 



For the case where the curve is a circle, equations flB.2[) - flB.4[) can be much 
simphfied. Denoting a = — f) as in section [3l the foUowing identities hold: 

, 2 - 

K = ^ = — sin a, d = —2tv sin ad 

d r = d • = cos(a), d • n = —d - n' = — sin(a) 
n ' n' = r ' r' = cos(2a), r - n' = —r' • n sin(2a) 

Substituting these in flB.2[) flB.4l) . and expressing the trigonometric coefficients in 
terms of c = cos a, the correlations become: 

Co = (l-c2)GPl + (l-2c2)Gwi 

d d^ 
C2 = (c^ - c2)GW + (12c^ - 12c2 + 

+ (28c^ - 32c2 + 5)G[2] 1 ^ (8^4 _ ;l2c2 + 4)G'[i] 1 (B.5) 

or, in terms of the actual derivatives G*^"^: 

Co = (1 - ^)G" + (-c')G'^ 

Ci = (c^ - c)G(3) + (3c3 - 2c)G"\ + (c^ - 2c)G'4 

a d'^ 

+ (7c4 - llc^ + 2)G"4 + (c^ - c2 + 2)G'-]- 
which is equivalent to fl31aj) fl31cj) . 
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Appendix C. Using asymptotic series to expand the integral transforms 

In section [321 and other sections, we discuss the number variance and form factor of some 
distributions. These statistics can be expressed as an integral involving the normalized 
correlation ( fl26l) and fl29l) ). for which we have an explicit expression. The expression 
for the correlation function (e.g. fl2Tl) or fl34l) ) is usually very complicated, which makes 
the resulting integral hard to solve directly. However, given the asymptotic expansion 
of the normalized correlation, we can calculate the asymptotic expansion of the integral 
up to a small number of constants, which can be calculated numerically. The outline of 
this method is described here. 

When calculating the number variance for large L, we encounter integrals of the 

form 

If{L)= r fix) Ax (C.l) 



JO 

where f{x) = x^TZ{x)^ and TZ{x) is the normalized correlation, for which the asymptotic 
expansion is known. We split the asymptotic expansion of / into two parts: 

no oo 

fix) r^^anix) + anix), 

n=l 72=720 + 1 

where no is the number of terms such that J^a^(x)dx diverges for every finite L 
(no may be 0, but we will assume it is finite. We will also assume that for n > Uq the 
integrals converge uniformly with some finite lower bound Lq). 
Defining 

720 

Dix) = J2 <^nix), Cix) = fix) - Dix), (C.2) 

n=l 

we find that C{x) generates an integral with known asymptotic expansion 

oo poo 

C(x)da;- Yl / anix)dx = - ^ 6„(L) (C.3) 

72=720 + 1 ^ 72=720 + 1 

We now wish to add and subtract D from the integrand of fIC.ll) . However, the 
interval [0, L] might contain points which would cause the integral of D to diverge. For 
such cases, we split the interval at some point Xo, which is larger than these problematic 
points, and do the subtraction only in the second half [xq^L]. We will assume that 
Xq and L are both larger than Lq. In practice, for the functions used in this paper, 
there is either no points of divergence (and Xq is chosen to be 0), or a single point of 
divergence at 0, in which case we choose Xq = 1 (which happens to simplify the resulting 
expression). Using flC.2[) . the integral becomes 



rxQ pL pL 

IfiL) = / fix)dx + / Cix)dx + / Dix)dx 
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The third term in this expression consists of the "diverging" terms of /, however 
the integral here is finite and its asymptotic expansion is easily obtained from flC.2[) : 

/ D{x)dx an{x)dx = ^ 6^(L) (C.4) 

To handle the remaining terms, we add and subtract flC.31) 

rxQ poo 

If{L) = / f{x)dx+ / C{x)dx 

Jo Jxo 

pL poo 

+ / D{x)dx - / C{x)dx 

J xq J L 

In this expression, the first two terms are finite values independent of L, and the 
last two have known asymptotic expansion. Using flC.31) and flC.4[) . we get the final 
result: 

oo 

If{L)^Mf + J2UL) (C.5) 

n=l 

where 

rxo 



PXQ POO 

Mf= f{x)dx+ / [f{x)-D{x)]dx 

Jo Jxo 
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